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Abstract 

Gapped excitation spectra of Andreev states are studied in one- and two-dimensional 
(ID and 2D) normal systems in superconducting contacts subject to a parallel 
magnetic field. In the ballistic regime, a specific interplay between magnetic field 
spin splitting and the effect of a screening supercurrent is found to preserve time- 
reversal symmetry for certain groups of Andreev states remaining gapped despite 
the presense of the magnetic field. In ID wires such states can lead to a fractional 
thermal magnetoconductance equal to half of the thermal conductance quantum. In 
2D systems the thermal magnetoconductance is also predicted to remain suppressed 
well below the normal-state value in a wide range of magnetic fields. 



1 Introduction 



Recently, there has been extensive experimental and theoretical work aimed 
at understanding various microscopic manifestations of the mesoscopic-scale 
superconducting proximity effect in normal metal (N)- superconductor (S) 
structures [1,2,3,4,5,6,7,8,9,10,11,12,13,14]. On the normal side of such sys- 
tems the superconducting correlations are maintained in the course of Andreev 
scattering [15] during which a particle incident from the N region with energy 
below the superconducting gap energy A coherently evolves into a Fermi sea 
hole with the opposite spin that retraces the particle trajectory (time-reversed 
path) back to the normal bulk. Due to elastic disorder in the N conductor or its 
restricted geometry, a large number of consecutive phase-preserving Andreev 
reflections can occur at the NS interface making it effectively transparent to 



pair current. This gives rise to a veriety of phase-coherent phenomena such as, 
e.g. the zero-bias enhancement of the electron conductance [1,2,3,4,5,11,16] 
above the value predicted by the theory of a single-event scattering at an 
NS boundary [17], the finite bias and magnetic field anomalies of the phase- 
coherent conductance [11,18,19,20], the formation of a superconducting mini- 
gap in the quasiparticle density of states in normal systems [21,22,23,24], An- 
dreev edge states [25,26,27] and bilhards [28,29,30,31,32]. 

While phase- coherent charge transport has been receiving considerable atten- 
tion, heat conduction properties of mesoscopic proximity structures have been 
explored to a much lesser extent. A few theoretical papers have dealt with the 
thermal conductance of Josephson junctions [33,34,35], Andreev barriers and 
interferometers [36,37,38], quantum wires with proximity-induced supercon- 
ductivity [39]. Only recently, there has been a breakthrough in experiments 
on thermoelectric properties of small metallic NS hybrids [40,41]. 

The purpose of this paper is to report a theoretical study of anomalous mag- 
netic field behaviour of heat transport in low-dimensional proximity struc- 
tures. The choice of low-dimensional systems is motivated by a progress in 
fabrication of superconducting contacts to high-mobility semiconductor quan- 
tum wells (see e.g. Refs. [42,43]). A number of experimental observations have 
suggested that due to the proximity effect such systems acquire properties 
of "clean" superconductors where the electron mean free path I is large com- 
pared to the induced superconducting coherence length .^^v [2,9,14,42,43]. More 
specifically, in planar semiconductor-superconductor junctions the proximity 
effect can be described in terms of Andreev bound states formed between the 
NS boundary and the back wall of the quantum well [2,21]. Each Andreev state 
is a mixed particle-hole excitation whose spectrum, according to the theory 
of Ref. [21], should have a superconducting minigap Eg smaller than A due 
to a residual interfacial barrier. The existence of the minigap can for instance 
provide a natural explanation for the unusually strong enhancement of the 
electric conductance reported in Refs. [2,30]. 

The ballistic character of electron motion in semiconductor quantum wells 
needs to be taken into account when studying a magnetic field influence on 
the proximity effect. As proposed in Ref. [39], a magnetic fleld B parallel 
to the plane of the quantum well [Fig. 1(a)] affects the Andreev states via 
a screening supercurrent induced in the superconductor. In this case a finite 
Cooper pair momentum 2P s{B) at the NS boundary violates the time- reversal 
conjugation of particles and holes which is accompanied by a Galilean energy 
shift pPs{B)/mN of the Andreev states {rriN and p are the electron mass and 
momentum in the plane of the quantum well). Since Ps{B) oc B, the energy 
shift leads to the gapless excitation spectrum and ultimately to the metallic 
behaviour of the thermal conductance. 
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The above mechanism of the magnetic field influence on the Andreev states ne- 
glects the energy of the Zeeman splitting ag^sB (where a = ±1/2, g and fiB 
are respectively the electron spin, g- factor and Bohr magneton). On the other 
hand, in InAs-based heterostructures which are commonly used for contacts 
with superconductors the g-f actor can be as large as 10-13 [42], and therefore 
the magnetic spin splitting may play an important role. At first glance, the 
combined influence of the screening supercurrent and the spin splitting on the 
superconducting proximity should anyway be a destructive one because both 
of them break time-reversal symmetry. However, as shown below the inter- 
play between the Zeeman energy ag/isB and the Galilean shift pPs{B)/'mr^, 
which are both linear in B, can lead to the existence of Andreev states with 
the minigap independent of the magnetic fleld. Unlike the case of spinless elec- 
trons [39], the thermal conductance remains anomalously small compared to 
that of a normal state even at relatively strong flelds, up to the critical flelds 
of the superconductor. 

The predicted behaviour of the thermal magnetotransport is also rather dif- 
ferent from the observed in diffusive superconductors [44,45,46] where the 
constructive interplay of the supercurrent and Zeeman effects is obstructed 
by strong momentum scattering due to which the linear Galilean term aver- 
ages out [47]. Instead, the orbital magnetic fleld influence is described by an 
isotropic in momentum space depairing energy which is of higher order in B 
(see e.g. Refs. [46,48,8]). Thus, the mechanism of the magnetic fleld influence 
on the Andreev states and their thermal conductance discussed in this paper 
is unique to ballistic proximity structures. 



2 Proximity effect on the excitation spectrum of ballistic electrons 

We first discuss the magnetic field infiuence on the proximity effect in a quasi- 
one-dimensional electron system (QIDES) coupled in parallel to a supercon- 
ducting film via a barrier of low transparency r as shown in Fig. 1(a). The 
generalization for two-dimensional normal systems will be given in the next 
section. The film thickness d is assumed much smaller than both the super- 
conducting coherence length and the London penetration depth in a parallel 
magnetic field B = [0; 0; 5]. 

As has been mentioned in the Introduction, in the geometry of Fig. 1(a) the 
subgap states in the quantum well are formed in the course of multiple Andreev 
refiections which mix particles and holes producing a minigap in the excitation 
spectrum [21]. For a weakly coupled quantum well, the minigap energy Eg ^ 
{vn/vs)tEo depends on the energy of the lowest occupied subband Eq, the 
interfacial transparency r -C 1 and the ratio of the Fermi velocities in the 
normal (vn) and superconducting {vs) systems. In a narrow QID channel the 
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Fig. 1. (a) Cross-sectional view of a superconductor (S)-normal system (N) junction. 
The N system can be either a QID channel along the y-axis or a two-dimensional 
electron system located in the y, z plane. The particle (p) and hole (h) momenta 
in the N system are both shifted by P5 = {e/c)Bd/2 in order to match the Cooper 
pair momentum at the surface of the S film, (b) Model for studying ballistic 
thermal transport through a proximity-affected normal channel (N) connecting two 
reservoirs with temteratures T and T + AT {AT < T). 

motion of the mixed particle-hole excitations along the channel (y-direction in 
Fig. 1(a)) can be desribed by two coupled equations for the annihilation and 
creation operators, ipapit) and xpl.a-pi't)'- 



ihdt 



2m N 



+ agiiBB + E, 



N 



-inot agfiB^ + -c/jv 



2mN 



^^,{t)=E,tar^l,_^{t), (1) 
^l,_p(t) = E,ta^"^'ij„,p{t), (2) 



Ps = {e/c)Bd/2, 



(3) 



where Egia2^' plays the role of the effective singlet pairing energy (cr2 is the 
Pauli matrix), Ej^ stands for the Fermi energy in the QIES, and p = is 
the quasiparticle momentum. The magnetic field influence on the quasiparticle 
spin is taken into account by the Zeeman term agfisB, whereas the orbital 
effect is described by the shift Ps of both particle and hole momenta due to 
the screening super cur rent induced at the surface of the superconductor [see 
also Fig. 1(a)]. As the thickness of the normal channel is considered neghgible 
compared to the superconductor thickness d, the shift of the electron and 
hole momenta in the QIDES can be taken equal to the surface Cooper pair 
momentum given per electron by Eq. (3) (e > is the absolute value of the 
electron charge) . Ps is proportional to the half-thickness of the superconductor 
reflecting the fact that the field fully penetrates the film and generates an 
antisymmetric (linear) distribution of the supercurrent density with respect 
to its middle plane. 



One should note that the description of the proximity-induced correlations 
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Fig. 2. Schematic view of the excitation spectrum (5) at finite magnetic fields 

{B = Bg): spin- up quasiparticles near p = pn (a) and spin-down ones near p = —pN 
(b) remain unaffected by the magnetic field due to the cancelation of the Zeeman 
and supercurrent effects at g = 2kNd. Dashed curves correspond to B = 0. 

by coupled (superconductor- like) equations of motion (1) and (2) has to be 
reconciled with the fact that in the normal system there is no intrinsic super- 
conducting pairing. As shown in Appendix A, the superconducting coupling 
in Eqs. (1) and (2) is induced through the boundary conditions at the NS 
interface, for which the thickness of the quantum well must be of the order of 
the Fermi wavelength. 

The solution of equations (1) and (2) is given by the Bogolubov transformation 
of the form 



V'ap(^) = Upbapex.p{-it€+p/h) + ia2' °'vpbi^_pex.p{-it€^p/h), 



(4) 



1 + 



vn{\p\ - Pn) 



K{\p\ - pnT + E^]y^ 



:4 



u. 



where hap and b^-a-p are Bogolubov's quasiparticle operators, Pn is the Fermi 
momentum and the excitation spectrum e^^ can be represented as 



^ap = ^NPssgnp - agfxsB ± [v%{\p\ - PnY + ^J]^^^ = 

= {k^dsgnp - ag)fXBB ± [vU\p\ - Pn? + E^YI^, (5) 

with kj^ being the Fermi wave-number. We neglect the term quadratic in Pg 
assuming Pg/2mN ^ E^. 

According to Eq. (5) the magnetic field shifts the quasiparticle dispersion 
curves at the points p = -^Pn with respect to the Fermi level. At a certain 
field Bg when Eg = {g/2 + kNd)iiBBg there appear first gapless excitations 
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near the Fermi points: spin-up ones at p — —pN and spin-down ones ai p = p^ 
as sketched in Fig. 2. This field can be expressed in terms of the fiux quantum 
$0) proximity- induced coherence length = fiVN/2Eg and superconductor 
thickness d as follows: 



Alternatively, it can be related to the critical fields of the superconducting 
film as: 



A E 

^ ^"''"U^ + g) ^ ^'''"-Ak^dii + gy 

where Bgrb ~ $o/ Arf is the critical field due to the orbital effect of the magnetic 
field (A is the London penetration depth) and Bgpin ~ A/fJ'B is the critical 
field of the paramagnetic Hmit [49]. According to Eqs. (7) for junctions with 
well separated gap energies -E^ -C A and kf^d g > 1 the field Bg is much 
smaller than both critical fields. 

It follows from Eq. (5) that for the special case where g = 2kNd the Zeeman 
and supercurrent terms cancel each other at one of the Fermi points and the 
minigap remains at the Fermi level [see Fig. 2]. For these excitations time- 
reversal symmetry is exactly preserved despite the presence of the magnetic 
field, li g 2kj^d, the excitations near both Fermi points eventually become 
gapless as the magnetic field increases. However in a finite range of fields given 

by 

Bg \^-g\ 



the minigap still exists at one of the Fermi points, which accounts for the 
anomalous behaviour of the thermal conductance of low-dimensional proximity 
systems discussed in the next section. 



3 Thermal conductance of low- dimensional proximity structures 



Since the pioneering works [44,50] studies of heat transport have been playing 
an important role in understanding spectral properties of low-energy excita- 
tions in superconductors (see e.g. Refs. [45,46]). In what follows we discuss 
how the unusual magnetic field behaviour of the quasiparticle energies in low- 
dimensional proximity superconductors refiect on their thermal conductance 
characteristics. 
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3.1 QID channel 



One more advantage of describing the proximity effect in the QIDES by the 
superconductor-hke operator equations (1) and (2) is that we can use the 
well known procedure of Ref. [51] to calculate the heat current jq. For a long 
ballistic QID electron channel between two reservoirs with temperatures T 
and T + AT (AT -C T) [see Fig. 1(b)], the formal derivation of jg was given 
in Ref. [39]. To discuss the combined influence of the Zeeman spin splitting 
and the supercurrent on thermal transport in the channel we can start with 
the following expression for jq [39] : 



P<PN 



P>PN 



J dpe^pdpe^pn<{e^p) + J dp e-^pdpelpn>{elp) 



P<-PN 



P>-PN 



(9) 



It is related to the excitation energies e^^ and the group velocities 9^6+^ of 
the "+" branch of the quasiparticle spectrum (5). We have taken into account 
the contribution of the "-" branch by exploiting the symmetry relation e~p — 
—^-a-p between the two spectrum branches. 

The four terms in Eq. (9) correspond to the two rightmoving {dp€^p > 0) and 
two leftmoving (i9pe+p < 0) quasiparticle modes of the excitation spectrum 
shown in Fig. 2. The distribution functions of the rightmovers and the left- 
movers are assumed to be set by the reservoirs, as n>(e^p) = n(e+p,T + AT) 
and n<(e^p) = n{e'^p, T), respectively, with n(e^p, T) being the Fermi function. 
The current (9) can now be expressed in terms of the differences n^{e'^p) — 
n<(ej^p), and finally the thermal conductance can be introduced as the pro- 
portionality coefficient between the heat current and the temperature drop, 
Jq = i^i{B,T)AT, where for k,i{B,T) we have 



TT 



/ 



2kBT 



cosh^ X 



+ 



oo 

/ 



X- 



'dx 



Eg+(ki^d—ag)ngB 
2kBT 



cosh X 



(10) 



Here Kin{T) = irkj^T/Sh is the thermal conductance of a normal channel [52] 
{kB is the Boltzmann constant). It is convenient to rewrite equation (10) in 
terms of the dimensionless magnetic field b = B/Bg, temperature t = ksT/Eg 
and the ratio of the Zeeman and supercurrent energies g [see equation (6)] as: 
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Fig. 3. Magnetic field dependence of the normalized thermal conductance for (a) 
ksT/Eg = 0.05 and (b) ksT/Eg = 0.5. 



iii{b,t,g) ^ _3_ 
Kin 



l-b 
. 2t 



cost? X 



l_(g-l)6/(g+l) 

2t 



/ 



l-(l-g)6/(l+g) 
2t 



x'^dx 
cosh^ x 



cosh^ X 



i+b 
It 



cosh^ x 



(11) 



In the absence of the magnetic field {b = 0) the conductance (11) is ex- 
ponentially small at temperatures below the minigap {t < 1). As shown in 
Fig. 3(a) low-temperature thermal transport can be stimulated by applying 
a magnetic field above the threshold (6) corresponding to the appearance of 
gapless excitations. However the character of the magnetic field behaviour of 
the thermal conductance crucially depends on the ratio of the Zeeman and 
supercurrent energies g. For g = 1 the conductance saturates at half of the 
normal metallic value Kin/^ in contrast with both cases of small and large g 
where it fully recovers the metallic behaviour. This difi'erence results from the 
compensation of the Zeeman and supercurrent effects for half of the excita- 
tions which remain gapped for ^ = 1 despite the presence of the magnetic field 
[see Fig. 2]. Figure 3(a) also shows that any small deviation from g = 1 even- 
tually drives the system to the metallic regime. In this case the anomalous 
" half- metallic" regime is present in a finite range of intermediate magnetic 
fields given by equation (8). If g is very close to 1, the upper hmit of this 
range B — Bg{l + g)/\l — g\ can approach the lowest of the critical fields of 
the superconducting film, Borb or Bspm- This means that despite small values 
of Eg compared to the superconducting gap energy A the proximity-induced 
thermal transport anomalies in one-dimensional systems can persist untill the 
order parameter in the superconductor is destroyed by the magnetic field. 
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Fig. 4. Normalized thermal conductance vs temperature for 6 = 5. 

Figure 3(b) demonstrates that the anomalous magnetic field behaviour of the 
thermal conductance can still be well distinguished at temperatures of order of 
Eg/ ks- This can be directly seen from comparing the temperature dependences 
for g = 1 and ^ = (or ^ > 1) as shown in Fig. 4. 

The parameter g in Eq. (6) characterizing the relative strength of the Zeeman 
and supercurrent effects involves the superconductor thickness rf, electron g- 
factor and the Fermi wave- number kj^. If the QIDES is formed in a semicon- 
ductor heterostructure, the latter two can be controlled by an external gate 
voltage. This in principle gives one more way of manipulating thermal trans- 
port in low-dimensional proximity structures. To investigate such a possibility, 
in Fig. 5 we present the dependence of the normalized conductance on the pa- 
rameter g. The surviving proximity effect manifests itself as a minimum of the 
ratio Hi{g)/ HiN around the critical point g — 1- This minimum reaches 1/2 at 
low temperatures [Fig. 5(a)] and remains well pronounced as the temperature 
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Fig. 5. Normalized thermal conductance vs ratio of the Zeeman and supercurrent 
energies g = g/2kNd for (a) ksT/Eg = 0.05 and (b) ksT/Eg = 0.3. Such a depen- 
dence can result from a gate voltage variation of the electron ^-factor or the carrier 
concentration (i.e. k^)- 
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approaches Eg/ks [Fig. 5(b)]. 



3.2 2D electron system 



Here we extend our approach to planar structures combining a superconductor 
and a two-dimensional electron system (2DES). The 2DES is chosen to be 
located in the y, z plane [Fig. 1] and has a form of a strip of width W (in the 
z-direction) with a large number of the electron channels: 'iW^pf^jh ^ 1. The 
supercurrent flows along the strip (in the y-direction) . The straightforward 
generalization of the excitation spectrum (5) for a 2D isotropic Fermi surface 
is 



^apoe = {kNd cose - ag)nBB ± [<(po -pjv)' + E^]^ , (12) 

where po is the absolute value of the quasiparticle momentum (po ^ Pn), is 
the angle between the quasiparticle momentum p = [po cos d\ po sin 9] and that 
of a Cooper pair P5 = [Ps; 0] in the 2DES plane. Note that in the direction 
of the supercurrent flow (6 = 0) and in the opposite direction {6 = vr) the 
quasiparticle dispersion (12) looks like that sketched in Fig. 2 near the Fermi 
points Po = pn and po = —pN, respectively. 

We assume that the temperature drop is created in the direction of the super- 
current flow, i.e. along the strip [Fig. 1(b)]. In this case the contribution of one 
electron channel to the net heat current is given by the ID ballistic formula 
(9). In the polar coordinates po, 9 the net heat current can be written as 



tt/2 



PO<Pl\ 



h? 



Po>Pn 

2pnW ^ J J dpoe^p^()n>{€:^p^())dp^€^\p,0+ 
+ / dpo e^p^g n< (e+ ^e) dp^ e+ 



+ 



PO>PN 



(13) 



where the integrals over the angles < 9 < n/2 and n/2 < 9 < n take 
into account the contributions of the channels with positive and negative Py, 
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respectively; dpyS'^lp^e is the quasiparticle velocity in the direction of the tem- 
perature drop which has to be expressed in terms of the polar variables po, 9. 
Note that with the assumed accuracy \Ps\ ^ Pn the supercurrent term in 
Eq. (12) does not affect the position of the minimum {pq ^ p^) of the quasi- 
particle energy e'^p^^. Calculating with the same accuracy the quasiparticle 
velocity, one finds dp^e+lp^e ^ dp^e^p^Q cos 9 . 

As in the ID case, to model ballistic non-equilibrium transport between two 
reservoirs with temperatures T+AT and T, in Eq. (13) we introduce different 
equilibrium distribution functions n^^ie'^p) — n(e^p, T -|- AT) and n^{e'^p) — 
n(e^p,T) for the quasiparticles with positive {dp^e'^lpf^e > 0) and negative 
i^py^alpod < 0) velocities, respectively. Now the integrals over the momentum 
Po can be transformed into the energy integration as follows: 



Jq 



2pnW 



7r/2 



cos 9d9 



de e \n^(e) — n^(e) 



Eg + {k^dcos 6—ag)iXBB 



2pnW 



h? 



71 LXJ 

^ j COS 9d9 j dee [n>(e) — n<(e)] 

7r/2 Eg+{ki^dcosO—ag)nBB 



(14) 



and the conductance K2{B, T) of the 2DES is introduced as Jq — K2{B, T)AT, 
where 



K2{B,T) 3 

K2n{T) tt 



n/2 



= ^ / COS 9d9 X 



2kDT 



cosh^ X 



+ 



oo 

/ 



x^dx 
cosh^ X 



After exchanging the integrals over 9 and x, the angle integral can be easily 
calculated: 



K2{b,t,g) 3 



E 



l-2agt/(g+l) 



2x'^dx 
cosh^ X 



(15) 



l-2agb/{g + l) 
2t 



+ 



I 



l-(2ag + l)6/(g+l) 
2t 



x^dx / _ [{2tx-l){g + l) + 2agbf \ ^ 
cosh^x V ) 
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Fig. 6. Normalized thermal conductance of a 2DES vs magnetic field for (a) 
kBT/Eg = 0.05 and (b) ksT/Eg = 0.5. 



(2ag-l)b/(g+l) 

J _^(h^l^_[{2tx-l){g+l) + 2agbf 



cosh^x V b'^ 

l-2ag6/(g+l) ^ 



2t 



where we again use the dimensionless magnetic field b, temperature t and the 
parameter g characterizing the relative strength of the Zeeman and super- 
current effects; k.2n{T) = {2pi^W/h)niN{T) is the thermal conductance of a 
normal 2D stripe. In principle, the 2D geometry allows for an arbitrary rela- 
tive orientation of the heat current and the supercurrent, and in this case the 
thermal conductance will depend on the angle between them. Our geometry 
where this angle is zero has an advantage of being equally suitable for studying 
heat transport in both QID and 2D systems. 

Figure 6(a) shows the magnetic field dependence of the normalized conduc- 
tance (15) for different values of g. Like in the ID case for small and large 
values of g the conductance quickly approaches its normal value when h exceeds 
1. However for g = 1 this increase is much slower due to the compensation of 
the Zeeman and supercurrent effects for spin-up quasiparticles moving along 
the supercurrent fiow {9 = 0) and for spin-down ones moving in the opposite 
direction {9 — ir). There is no saturation in this case because the majority of 
the excitations (i.e. with the intermediate angles 9 ^ 0,7t) are still affected 
by the magnetic field: they all become gradually gapless as the field increases. 
The nonmonotonic dependence of the normalized conductance on the param- 
eter g shown in Fig. 7 also suggests that for ^ pa 1 the proximity effect in the 
2DES survives in higher magnetic fields. 
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Fig. 7. Normalized thermal conchictancc of a 2DES vs parameter g = g/ik^d for 
(a) kBT/Eg = 0.05 and (b) ksT/Eg = 0.4. 

4 Conclusion 



An intermediate regime, between superconducting and metallic, of heat trans- 
port has been predicted in baUistic low- dimensional proximity superconduc- 
tors subject to a parallel magnetic field. It is due to the combined influence 
of the screening supercurrent and the Zeeman spin splitting on the proximity- 
induced gap (minigap) in the quasiparticle states. The interplay of the two 
effects is found to preserve the minigap at the Fermi energy for certain groups 
of the quasiparticles despite the presence of the magnetic field. In this regime 
the low-temperature thermal conductance is about two times smaller than in 
the normal state and nearly independent of the magnetic fleld. For compari- 
son, in conventional superconductors with a supercurrent [46,47,8] or the spin 
splitting [46,49] acting separately, no gapped excitations remain at the Fermi 
energy as the field increases, and the thermal conductivity fully recovers the 
metallic behaviour. 

To observe the anomalous behaviour of the thermal conductance the param- 
eter g/2kNd controlling the relative strength of the Zeeman and supercurrent 
effects must be close to 1. This requirement can be met in nanostructures 
combining large g-factor semiconductor quantum wells (e.g. InAs-based) and 
thin superconducting films with typical thicknesses d ~ A;^^ ~ 10 nm. It is 
important to emphasize that the magnetic field Bg [Eq. (7)] characterizing 
both supercurrent shift and Zeeman splitting of the Andreev states is small 
compared to the critical fields of the superconducting film (arising from the 
orbital or paramagnetic effects). For this reason, we did not consider the mag- 
netic field dependence of the order parameter in the superconductor and the 
effective pairing energy Eg in the normal system. This question as well as the 
role of elastic scattering could be a subject of further studies. 
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A Effective pairing energy description of the superconducting prox- 
imity effect in low-dimensional systems 

Below we present the microscopic derivation of the equations of motion (1) and 
(2) and discuss their applicability. The approach used here is close in spirit to 
the tunneling self-energy description of the superconducting proximity effect 
proposed in Ref . [53] . However we will not employ the tunnehng Hamiltonian 
formalism [54] whose validity for high-order tunneling processes was a subject 
of debates [55,56,57]. Instead, to generate the tunneling self-energies we will 
use directly the equations of motion in the superconductor in the position 
representation. It turns out that in this approach there is no need to resort 
to the model of a homogeneous pairing potential in the superconductor. In 
this sense, a more general treatment of the proximity effect can be given and 
the results can then be compared with those obtained for the homogeneous 
pairing [21]. 

We assume a rectangular potential barrier with flat walls occupying the region 
—a < X < which separates a superconductor {x < —a) and a normal system 
{x > 0) [see Fig. 1(a) where the barrier is shown as a thin black rectangle]. 
The Fermi energies Es^n, Fermi momenta ps,N and the electron masses ms^N 
in the S and N systems are considered different. For simplicity, the derivation 
will be given for zero magnetic field. 

The superconductor is desribed by the mean-field equations of motion for the 
creation ipa-p{x) and annihilation ipapi^) operators [58,59]: 

[ihdt - es{Px, p) + Es] ipap{x, t) + Aa'a{x)ipi'-p{x, t) = 0, (A.l) 
[-ihdt - esi-Px, -p) + Es] ^l-pix, t) + A*^,^{x)jlja'p{x, t) = 0, 

where Aa'ai^) is the pairing potential, es{px, p) is the one-particle operator of 
the kinetic energy with p = \py,Pz] being the parallel (conserved) component 
of the electron momentum and Px — —ihdx, and we assume a summation over 
the spin projection a'. In the normal system (initially three-dimensional) the 
equations of motion are 

[ihdt - £n{Px, p) - V{x) + En] ipap{x, t) = 0, (A.2) 
[-ihdt - £n{-Px, -P) - V{x) + En] '0l-p(a;, t) = 0, 
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where Sn{px,p) is the operator of the kinetic energy of a normal electron, 
V{x) is a confining potential. Inside a high enough barrier one can neglect 
the energy and momentum dependence of the electron penetration length and 
write the equation of motion as: [dl—q'^]Xap{x) = 0, where = h/ {2mBUy/'^ 
is the electron penetration length depending on the barrier height U and the 
electron effective mass rriB- We introduce a special notation Xap{x) for the 
electron operator inside the barrier. The continuity of the particle current 
imposes usual boundary conditions at the barrier walls: 



Xap(O) = ^ap(O), Xap(-O') = ^ap(-a), (A.3) 

dxXap{0) = —dxi/jap{0), dxXap{-a) = —dxi^api-a)- (A.4) 
rriN ms 



The solution inside the barrier satisfying boundary conditions (A.3) is 



. . sinhg(a; + a) , sinhga; , , 

Xap{x)^ — V'ap(O) - ^-r — Wap{-a). (A.5) 

smhga smhga 

Inserting it into the boundary conditions for the derivatives (A.4), we have: 



— 5x'0ap(O) - goV'ap(O) = -gt Vap(-a), (A.6) 

rriN 

—dxipap{-a) + qoipap{-a) = QtipapiO), (A.7) 
ms 

Qq — q/ tanhqa, qt — q/sinhqa. (A. 8) 

These equations serve now as effective boundary conditions for the equations 
in the superconductor and the normal system. For a thick barrier where qt — > 
0, the coupling between the "normal" and the "superconducting" operators 
vanishes, which is described by Eqs. (A.6) and (A.7) with zero right-hand sides. 
Using Eqs. (A.6) and (A.7), one can make sure that the current is continuous 
at the interface: 



Jap 



ih 
2m N 



5.^1p(0)V'ap(0) - h.C 



ih 



2mn 



2ms 
^ip(0)^ap(-a) 



dx'4^ip{.-o)'4^ap{-o) - h.C 



h.C 



In what follows we assume that the influence of the normal system on the 
superconductor is negligible and the proximity effect in the normal system can 
be studied without a feedback. It is normally the case if the superconductor is a 
bulk metallic system whereas the normal system is a degenerate semiconductor 
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where the electron concentration is much lower than in the metal. In terms of 
the Fermi momenta this assumption can be expressed as: pjy ^ Ps- Besides 
we restrict ourselves to the energies e much smaller than the gap energy A 
in the superconductor: |e| ^ A. In this case, we will derive closed boundary 
conditions for the normal electrons which will take into account the proximity- 
induced electron pairing. 

It is convenient to include the effective boundary condition at the S side of 
the barrier (A. 7) into the equations of motion in the superconductor (A.l) by 
adding appropriate delta-functional terms as follows: 



ihdt - esiPx: P) - Us(x) + Es] ^l^apix) + Ac,'a(a;)V^i,_p(a;) = 
= -|^<^(^ + a)V'ap(0), 
-ihdt - £si-Px, -p) - Us{x) + Es ^i-p{x) + A*^,^{x)^a'p{x) = 



(A.9) 



2m F 



S{x + a)jjji_p{0). 



where a singular potential 




reproduces the boundary condition (A. 7) with zero right-hand side which cor- 
responds to an isolated superconductor. The penetration of Andreev bound 
states in the superconductor at low energies is described by a particular solu- 
tion of Eqs. (A.9) generated by the right-hand sides containing the operators 

of the N system. It can be expressed in terms of the matrix Green function 
of Eqs. (A.9) whose matrix elements are constructed from the quasiparti- 
cle Gp" {x, t; x', t') and condensate (Gorkov) Fp" (x, t] x\ t') Green functions, 
namely: 



^'^i-p(x,t) ^ 



2mf 



X 



(A.ll) 



X 



{x, t- -a, t') (x, t- -a, t') \ ( ^„,p (0, t') 

{x, t- -a, t') Gif {x, t- -a, t') ) [ (0, t') 

for an isolated 



These Green functions can be found from Gorkov equations 
superconductor. [The star means complex conjugation] . 
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Inserting solution (A. 11) into the effective boundary conditions for ipap and 
■0a-p the N side of the barrier [Eq. (A. 6)], one finds 




(A.13) 



(A.12) 



For low-energy states (|e| << A) these equations are reduced to: 



TUB 




GV'ap (0) + F*i<7?">l'-p (0)1 , 



2mB 



(A.14) 



rriB 




CVLp (0) + FzarVa'p (0)1 , 



rriN 



2m,B 



where the Green functions are taken at the surface of the superconductor at 
zero energy and momentum: G = G{—a; — a)|e,p=o and F = F{—a; — a)|e,p=o- 
The dependence on the parallel momentum is ignored because for |p| < ^ 
Ps only electrons moving nearly perpendicular to the interface can tunnel from 
the superconductor into the normal system. 

Boundary conditions (A.14) take into account the conversion of a particle into 
a hole (and vice versa) due to Andreev reflection that occurs simultaneously 
with normal scattering. The advantage of using equations (A.14) is that the 
Andreev process is described by an explicit coupling between the ^ap(O) and 
■0-a-p(O) operators of the N system with the condensate Green function of 
the superconductor F (and also G) reduced to a constant. Another advantage 
is that from the boundary conditions for the operators (A.14) one can easily 
derive the boundary conditions for both Green functions and Bogolubov-de 
Gennes functions of the N system. 

To study the effect of superconducting proximity on the excitation spectrum 
of a clean N system one needs to take into account its flnite size determined 
by the confining potential V{x) in the equations of motion (A. 2). Due to 
the back wall of the quantum well particles and Andreev reflected holes are 
scattered back to the NS boundary and interfere with those outgoing from 
the interface, which generates (in a very subtle way) an energy gap in the 
quasiparticle spectrum [21]. However in a very narrow well (two-dimensional 
electron system) where the electron states are localized near the NS surface 
the proximity-induced correlations can be treated in a more robust way. In 
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this case the couphng between ii'„p(0) and '?/'-o-p(0) in the boundary condi- 
tions (A. 14) represents an effective in-plane pairing for the two-dimensional 
electrons. Indeed, combining equations of motion (A. 2) and boimdary condi- 
tions (A. 14) and neglecting all the terms not involving the condensate Green 
function, one can write: 



[ihdt - SNiPx, p) -V{x) + En] ipap{x, t) = 

(A.15) 

[-ihdt - £n{-Px, -p) - V{x) + En] ^la-p{x, t) = 
= <5(^)(|£) Fi<jr"Va'p(^,i)- 

Taking into account that for a weakly coupled 2DES the transverse wave func- 
tion ip{x) is almost unaffected by the tunneling, one can multiply Eqs. (A.15) 
by ip{x) and integrate over x (i.e. over the thickness of the 2DES). This leads 
to the equations of motion (1) and (2) for the operators of the lowest occupied 
subband with the effective pairing energy Eg given by 



To estimate Eg one can use the condensate Green function of a supercon- 
ductor with a homogeneous pairing potential A at zero energy and parallel 
momentum: F ^ J2p^ ^ / {^'^ + i^siPx — PsY)- The integration over the 
perpendicular momentum gives F ^ l/hvs- The boundary value V'(O) of 
the transverse function can be estimated using boundary condition (p{0) ~ 
{tub / mN(l)dx^{S^) [see Eq. (A. 6)], where in the right-hand side one can use 
the "hard wall" wave function </?(x) = {2/ Ly^'^ sinnx/ L, which gives </?(0) 
{niB / 'mNq){2/ LY^'^{'k / L) [L is the thickness of the quantum well]. Thus, the 
effective pairing energy is 



Eg = ^ — ~ — T^o, A.17 

niNJ-iVs smh qa 2m n^^ vs 



with r = 1/ sinh^ qa <^ 1 and the energy of the lowest occupied subband 
given by Eq = h'^Ti'^ /2mNL'^- Equation (A.17) is equivalent to that obtained 
in Ref. [21] for a strong delta-shaped barrier. 
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